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A SIMULATOR INVESTIGATION OF ENGINE FAILURE COMPENSATION FOR 
POWERED-LIFT STOL AIRCRAFT 
Albert W. Nieuwenhuuse 
and 

James A. Franklin 

Ames Research Center 
Moffett Field. California 94035 

ABSTRACT 


A piloted simulator investigation of various engine failure compensation concepts for powered- 
lift STOL airaaft was carried out at the Ames Research Center. The purpose of this investigation 
was to determine the influence of engine failure compensation on recovery from an engine failure 
during the landing approach and on the precision of the STOL landing. The various concepts 
included ( 1 ) cockpit warning lights to cue the pilot of an engine failure, (2) pro«rammed thrust and 
roll trim compensation, (3) thrust command and (4) flight-path stabilization. The aircraft simulated 
was a I 50 passenger four-engine, exte'nally blown flap civil STOL transport having a 90 psf wing 
loading and a .56 thrust to weight ratio. Results of the simulation indicate that the combination of 
thrust command and flight-path stabilization offered the best engine-out landing performance in 
frbulence and did so over the entire range of altitudes for which engine failures occurred. 
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ON THE QUESTION OF PILOT DETERIORATION DURING LOW ALTITUDE FLIGHT 

P, Schulz* 


1. Introduction 


/ 99 » 


The ride qualities of an aircraft represent one of several 
multi-discipline prob]eras which are related to the development 
of a controller for the high velocity low altitude flight regime. 

The usual criteria for evaluating the ride qualities of an 
aircraft are the variances of accelerations at the location of 
the pilot seat. Since the factors which Influence the minimiza- 
tion of the aircraft response to gusts also Influence the flight 
properties of the system, we should also consider the following: 

1. Stability and control, and 

2. Instabilities Induced by the pilot, 

l.e., a gust aid system cannot be designed in a manner which is 
too stable, otherwise, the pilot will have difficulties in 
flying tl.e aircraft. 

However, these problems will not be discussed in this paper, 
which is only concerned with the calculation of the ride qualities 
during the project and definition phases. In ordei to be able to /lOO 
calculate the ride qualities in a satisfactory manner, the accur- 
acy of the following models is Important: 

*VFW-Fokker GmbH., Bremen. 

**Numbers in the margin indicate the pagination of the 
original foreign text. 
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1. the turbulence model for low altitudes,. 

2. the aircraft model, 

3. the aircraft mission model, 

4. the pilot model which is loaded by the accelerations. 

The first calculations can be made using a I'igid aircraft. 
Later on, when the elastic and aerodynamic properties are 
better known, it is possible to Include the elastic degrees of 
freedom. There are two possibilities for calculating the trans- 
fer function of the aircraft: 

1. the equations for longitudinal motion of the point 
mass aircraft, 

2. buildup of a relatively simple discrete aircraft model 
and use of the rigid degrees of freedom for calculating the 
transfer function. 

Both methods will be discussed because they give different 
results for aircraft with large sweepback ang?es. 

The tolerance limits of human beings to accelerations have 
not yet been determined even today, because not enough data have 
been collected. 

The pilot deterioration should be the general criterion for 
the ride qualities, and not the variance of the accelerations at 
the pilot seat. This requirement is a consequence of the fact 
that two different power sptctra can have the same variance 
(because of the in^ er-'At I :n j but could result in different degrees 
of pilot deteriorar ■i . Therefore, this lecture will primarily 
deal with this topic. 
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The sensitivity of an aircraft to atmospheric turbulence 
primarily depends on the ratio of the force produced by the gust 
and the mass of the aircraft, as well as the tendency of the 
aircraft to turn into the relative wind, In order to calculate 
the gust sensitivity at the center of gravity and at the pilot 
seat, we can use the equations for the longltu'^lnal motion as a 
first approximation [3l. The results of these calculations will 
be very useful for the following investigations: 

1. the influence of the angle of attack oscillation on 
the gust sensitivity at the center of gravity. The result is 
that there is a decrease in when ^ and are increased 

(Figure 1). 


2. the influence of the angle of attack oscillation and 
of the dynamic amplification factor on the gust sensitivity at 
the pilot seat (Figure 2). 

Since these various parameters for optimum ride qualities 
also Influence the flight property criteria, for example , 

proper selection of these quantities 
depends on both the ride qualities and the flying qualities. 

2. The Behavior of Humane under the Influence of 
Stochastic Disturbances 


Not much data is available which describes the fatigue of 
humans caused by stochastic accelerations. When these data are 
applied to high velocity, low altitude flight, one encounters 
the problem of whether and how the pilot can solve tactical and 
flight handling tasks If he is subjected to a stochastic accel- 
eration spectrum. All of the data on this subject apply for 
sinusoidal accelerations with a certain variance and frequency, 
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M« = a9 



Figure 1. Influence of angle of attack oscillation on the gust 
sensitivity at the center of gravity. 

From Ralph C. A. Harrah. Maneuverability and Gust Response 
Problems Associated with Low Altitude, High Speed Flight. 
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Figure 2. Influence of angle of attack oscillation on the gust 
sensitivity at the pilot seat location. 

From Ralph C. A. Harrah. Manueverabllity and Gust Response 
Problems Associated with Low Altitude, High Speed Plight. 



for example, see [1]. This will be used here as a basis for 
calculating the pilot deterioration due to gusts. This reference 
indicates that the acceleration within six octaves must be 
Integrated to obtain partial values of the variance, l.e., 

UK 

W; 

where^^(fco)= power spectrum of accelerations. Using these 
partial values of the variance of the accelerations, the per- 
missible loading times t are calculated from the load capacity 

^1 

curves (Figure 3). The reciprocal S. «= 1/t is then the pilot 

1 

deterioration. The deterioration from all octaves Is then summed 
and one finds the total deterioration due to the existing power 
spectrum of accelerations at the pilot seat location. 


The method given In this reference uses a type of average 
value of accelerations in this octave for calculating the 
-’eterloratlon . Since the permissible loading times depend very 
greatly on the variance and frequency of the accelerations, the 
pilot deterioration will change drastically when the octaves are 
displaced (Figure ^). Therefore, it would be desirable to have 
a method which does not use the average value of the accelerations 
but a method which averages the permissible loading times for the /103 
deterioration. This can be done by means of the following substi- 


tution 


U> » to- 


■ Z 


Jilo » Z ’ u>^ 2 Jly 
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In this formula, x la a dlmenslonlesa quantity, l-.f., the 
Integrand has the same dimensions as o*. 


Therefore, we make the following substitution 

- w , 2 ’‘ J 8*.2 

and we obtain the following for the variance of octave 1 

“ 1 


This expression has the advantage that 
function and has the same dimensions as o* 

part , 


Is a continuous 
Theref ore , 


^*part average of the function of the octave 1, 

because the length of the Integral Is 1. Therefore, the 
deterioration curve S*(a*, x) Is calculated from the continuous 
j and the permissible loati curves (Figure 3). The 


function 


earlier summation of partial deterioration values is then 
replaced by an integration of the deterioration curve. 

^ J d[y 

The disadvantage of the old method is that it averages 


This fact and the variation of the pei’mlsslble load curves /104 
are the reasons for the change In the results when the octaves 
are displaced. 
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The question of tie applicability of the permissible load 
curves to stochastic accelerations becomes rather unimportant by 
making a transfer to the continuous function from the partial 

values of variance. These permissible load curves are then es- 
tablished for sinusoidal accelerations with a certain frequency 
and amplitude. Discrete points of the curve i.e., the 
corresponding power spectrum of accelerations, can be looked upon 
as representations of sinusoidal functions with frequencies and 
variances corresponding to the selected points. Therefore, the 
reference mentioned above can also be applied for stochastic 
accelerations 

3j Pilot Deterioration Due to Stochastic Accelerations 

A pilot can withstand an acceleration spectrum with a 
certain variance for only a limited time t^, which is a function 
of the average value of accelerations and the frequency of the 
gusts. This function was taken from [1] mentioned above. The 
variance of the load multiple o at the pilot seat is given by 

09 ” 

0 

00 

is the power spectrum of the load multiple; (^^(u) is 
the power spectrum of the gusts; P(w) is the transfer function /10 3 
of the aircraft Avi^/w . 

During the mission, the transfer function of the aircraft 
changes because of the decrrase in aircraft weight. Other reasons 
for this are the changes in weather, velocity, and flight altitude 
during the flight. In order to include these changes in the 
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load multiple In the calculation, the Influence of any effective 
time Is explained by means of a damage accumulation [2]. 


The "damage rate", which varies In time, has been glv * by 

ii J— 

where t (t) = instantaneous permissible load time due to a (t). 

The Instantaneous damage Is then the ratio of the loading time 

--dt — and the permissible loadlnr time t (t). 

z 


In this way, the total damage over the mission becomes 


S-^CTJ- 

where Tj^^ = mission time. 




Since any mission differs from any other mission because of 
the ehange In the weather and because only the statistical distri- 
bution of the gusts is known. It is Impossible to calculate the 
actual damage during a certain mission. Therefore, it Is only 
possl'le to calculate the expectation value of pilot deterioration: 




S(<r») JlVKo-w) 
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where W(a ) = probability function and w(a ) = probability 

W W 

density function. 


It is defined in [^] as 





follows : 




11 


i 


where the parametei-s and are the ratio of flight time or 

flight path for nonstormy and stormy turbulence, and b^ are 

the jcale parameters for the Individual probability distributions 
of 0 ^ for the two types of turbulence. 

The power spectrum for the vertical turbulences is given by 
the Dry den or the Karmann spectrum. 

The Calculation of the Aircraft Transfer Function 


As already mentioned '<bove, there are two possibilities for 
calculating the transfer function of a rigid aircraft. 

1. The equations of motion for the point mass aircraft 

2. If a relative simple discrete model Is used, it is 
possible to employ the rigid body motions for calculating the 
transfer function. 

The main differences between tVie method are In the treatment /107 
of the unsteady aerodynamic effects and the sweepback angle. 

In order to calculate the transfer function of the rigid 
aircraft due to vertical gusts for the longitudinal motion, the 
equations of motion are taken from [5] 

is ~ ^ 

The unsteady aerodynamic effects are not contained in these 
equations. Since the deflections of Individual Wxng points 
cannot be determined in this model, it is only possible to 
consider the Kuessner function. Therefore, for the right side, 
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t 


we obtain 


t 

where \P - Kuessner function. 

The transfer function is obtained from the two equations 
given above by using the Laplace transformation. The pitch 
moment term In the moment equation can be split up Into compon- 
ents for the wing-fuselage and for the elevator, in order to 
Include the dead time between these two components. The Kuessner 
function of the elevator can be ignored because the chord length 
of the elevator is small compared with that of the wing and the 
aerodynamic response can the'^fore be approximated by a jump 
function. 

The transfer functions for the two degrees of freedom of the 
system a and ^ are calculated from the two equations of motion. 
With this it is possible to then calculate the load multiple at 
the center of gravity and the additional load multiple caused by 
rotation around the center of gravity. 



The transfer function for the total load multiple is then the 
sum of both 

U ’ W ^ w 
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The calculation of the transfer function of the elastic 
aircraft will be outlined here. Starting with the general matrix 
equation of an elastic equation^ It is possible to determine the 
displacements and rotations of the discrete points of the air- 
craft, These are then built up from the contributions of the 
individual eigen modes and of a time function, according to the 
Bernoulli separation theorem. When a transfer is made from the 
time region to the frequency; region, this then becomes the trans- 
fe:- function of the system. Because of the fact that the deflec- 
tions are known, it is possible to apply the Kuessner and the 
Wagner force under these unsteady aerodynamic conditions, In our 
calculations, we only considered the circulatory components 
caused by vertical motion of the profile and of the rotation 
around the C/iJ point in the Wagner force. 

The advantage of an aircraft built up of discrete points Is 
that it also becomes possible to consider the influence of the 
wing sweepback. In the case of a wing with a strong sweepback, 
a straight gust front will reach the Individual points of the 
wing at different times and, therefore, the response of the 
aircraft will be softer. Another advantage Is the possibility 
of being able to consider the Kuessner and the Wagner forces 
under unsteady aerodynamic conditions. 

These two facts are the reasons for the differences In the /109 
two curves given In Figure 5* This is also the reason for their 
general appearance. Even though one wishes to only calculate 
the rigid aircraft, the discrete model can be used because the 
rigid degrees of freedom and the eigen modes can be used for 
describing the aircraft and its motion. In addition, later on, 
the elastic eigen modes can be incorporated relatively easily 
into the computer program. 
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Summary 


The calculation of pilot deterioration for high velocity, 
low altitude missions is a multi-layered problem. The Improved 
method of considering the physiological data of humans represents 
only a small step towards the exact calculation of ride qualities. 
Because of the difficult situation of a pilot during a terrain- 
following mission, it is necessary to have a flight controller 
for controlling the aircraft. Therefore, a complete aircraft 
model for calculating the pilot deterioration should contain this 
controller . 

The load on the pilot caused by maneuver stresses during 
terrain following flights is very hlgi). Therefore, it is 
necessary to consider these accelerations in the calculation. 

However, we encounter the following difficulty: if one super- 

imposes the acceleration caused by vertical gusts with those of 
the terrain following flight, one then Indirectly is saying that 
there is no relationship between the gusts and the terrain shape. 

Even if we have a noise generator and a filter during a real 
time simulation, which produces the statistical gust distribution 
along the flight path, one has still not yet made any statement 
about the relationship between these two quantities. Another 
difficulty is the low frequency at which the maneuver loads of 
terrain following flight occur. The permissible load curves given 
in [1] extend to a frequency of down to f = 0.7 Hz. After about /IIP 
1 Hz, one encounters the region of seasickness, which was inten- 
tionally excluded in this report, because not enough is known and, 
in addition, the available data is contradictory. 

For the first calculations of pilot deterioration during the 
project phase, it is not necessary to consider the elastic degrees 
of freedom. Since the influence of the elastic oscillations on 
the pilot is considerable even for small aircraft, one should 
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include the elastic degrees of freedom as soon as possible In the 
calculations. The curves gi^-en in Figure 6 show the Influence 
of the elastic degrees of freedom. The slight difference between 
the two curves at large swe« pback angles was produced by the 
fact that the moments of inertia of the wing segments were too 
small . 

The method for calcula Ing the pilot deterioration is the 
same for sideways motion. Only the permissible load curves are 
somewhat different in this case. 
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Figure 1. Time History of STOL Landing 

Following an Engine Failure at 50 Feet 
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Figure 2. Trends of Engine-Out Touchdown Sink Rate with 
Altitude at Engine Failure 
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Figure 4. 
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Figure 6. Flap Assembly and &;gir® Pylon Detail 
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Configuration - ALL Engine? Operating 
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Figure 10. Longitudinal Response to a Step Column 
Input - SCAS Off and On 
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Figvure 12. Lateral-Directional Response to 
Wheel and Pedal Doublet Inputs - 
SCAS Off 
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Figure 15. Engine Thruat and RPM 
Relation ship a 

























3.0 rmd/uc 


O 



Figxire 19. Root Locus and Bade Plots of 
Closed-Loop Thnist Control 
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Flgiire 20. Transient Thrust Response to a 
Thrust Command Step Input - 
Closed Loop Thrust Control Off and On 
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Figure 22. Time History of a Landing Flare Performed 
with Aircraft Rotation - All Engines 
Operating 
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Figure 25. The Effects of Warning Lights 

on Ekigine-Out Landing Performance 







- Engine Failed at 86 Feet 






Figure 27 (concluded) 
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Figure 28 (concluded) 
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Figtffe 29. The Effect of Programmed Thruat 
and Roll Trim Compeneatlon on 
Engine-Out Landing Performance 
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Figure 29 (concluded) 
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Figrire 30. Longitudinal Response to a Step Column 

Input for Several Flight-Path Stabilization 
Feedback Gains 
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Figure 31, Tiae History of a Landing with Flight Path 
Stabilisation I^a3p®rcent per degree 
-ill Engines Operating 
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Figvire 32, Time History of a Landing with Flight Path 
Stabilization ^ = 5 percent per degree 

h- 

- All Engines Operating 
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Figure 33. Time History of a Landing with Flight Path 
Stabilization,!^ s/0 percent per degree 

h 

- All Engines Operating 
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Figure 34. Comparison of Landing Precision of Basic 

Aircraft with and without Flijht-Path Stabilizatic-i 
- All Engines Operating 
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Figure 35. Comparison of Touchdown Sink Rate 

Probability for Aircraft With and Without 
Flight-Path Stabilization 
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Figure 36 . Time History of an Engine-Out Landing 

- Flight Path Stabilization, s 3 percent per degree 

- Engine Failed at 61 Feet ~ 
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Figure 37. Time History of an Engine-Out Landing 


- Flight Path Stabilization, fc^B 5 percent per degree 


- Engine Failed at 86 Feet 
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FigLire 38. The Effect of Flight-Path Stabilization 
on Engine-Out Landing Performance 
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Figure 39. Time History < 

- Closed Loop 

- Engine Failed 
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Figure 40« The Effect of Cloaed-Loop Thrust 
Control and Combined Cloaed-Locp 
Thrust ControlAllght Path Stabilization 
on Englr.e-Out Landing Performance 
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Figure 41. Time History of an Engine-Out Landing 

- Closed-Loop Thrust Compensation and 

Flight Path Stabilization, 5 percent per degree 

- Engine Failed at 86 Feet ^ 
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Figure 42. Comparison of Trends of Engine-Out 

Touchdown Sink Rate with Altitude of 
Engine Failure for the Engine Failure 
Compensation Concepts Investigated 


